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When magnetic field is applied parallel to the ab plane in d^2_y2-wsLye superconductors, the transition of 
stable vortex lattice structure, spatial structure of local density of states, and specific heat oscillation by ro- 
tation of magnetic field orientation are investigated by quantitative calculations based on the selfconsistent 
Eilenberger theory We estimate how the vortex state changes depending on the relative angle between the 
node-direction of the superconducting gap and magnetic field orientation. To reproduce the sign-change of 
specific heat oscillation observed in CeCoIns, our study is done by including strong paramagnetic eff'ect. 
The quantitative theoretical calculations give decisive information to analyze the experimental data on the 
field-angle dependence, and establish the angle-resolved specific heat experiment as a spectroscopic means 
to identify the node-position of the superconducting gap. 

KEYWORDS: vortex state under parallel field, paramagnetic effect, specific heat oscillation, Eilenberger 
theory 



1. Introduction 

In unconventional superconductors such as heavy fermion 
superconductors, high-Jc cuprate superconductors, organic 
superconductors, the pairing mechanism may not be con- 
ventional electron-phonon interaction. The Cooper pairs are 
formed by exotic pairing mechanism coming from the strong 
electron-electron interaction, and possibly the pairing symme- 
try is p-, d-, or /-wave, other than full gap ^-wave. Thus, the 
pairing function has sign change on the Fermi surface, result- 
ing in nodes of the superconducting gap. To clarify the mech- 
anism of unconventional superconductors, it is important to 
identify the pairing symmetry of the Cooper pairs. 

Experimentally, existence of the node can be detected by 
the power-law dependence of physical quantities as a function 
of temperature T. For example, the specific heat C oc (oc 
T^) and nuclear spin relaxation rate oc (oc T^) for line 
(point) node of superconducting gap, instead of exponential 
T-dependence for full gap 5-- wave pairing. This comes from 
the power-law of the density of states (DOS) N(E) oc E (oc E^) 
for the line (point) nodes of the gap function. 

The position of the node on the Fermi surface can be identi- 
fied by the phase- sensitive experimental methods such as cor- 
ner junction^' and Andreev surface bound states."^^ For these 
experiments, we need delicate fabrication techniques to set 
up a good junction. Another method to detect the node po- 
sition is magnetic field orientation- sensitive experiments for 
physical quantities of bulk measurement. There, by rotat- 
ing magnetic field orientation, we try to identify the position 
of the nodes from the change of physical quantities, such as 
specific heat,^~^^^ and thermal conductivity, depending on 
the relative angle between magnetic field orientation and the 
node-direction. These quantities under magnetic fields reflect 
electronic excitation in the vortex states. For reliable analy- 
sis of the experimental data and establishing those methods 
as a spectroscopic means of the gap node position determi- 
nation, these electronic behaviors have to be confirmed by 
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Fig. 1. (Color online) Schematic configuration of magnetic field and 
<i_^2_-y2-wave superconducting gap on rippled cylindrical Fermi surface. We 
discuss the case when magnetic field is applied parallel to the ab plane. 



quantitative estimate of theoretical calculation in the vortex 
states. In this paper, we study the vortex states and the field- 
angle sensitive specific heat under magnetic fields parallel to 
the ab plane, related to the pairing gap in a heavy fermion su- 
perconductor CeCoIns. For CeCoIns, there were controversial 
discussions as for the pairing symmetry between J;c2_^2-wave 
and Jxy-wave.^^^ 

We consider the case when the magnetic field orientation 
is rotated within the aZ^-plane, and the Fermi-surface is quasi- 
two dimensional (Q2D), as shown in Fig. 1. In this ab plane 
magnetic field configuration, we calculate the vortex structure 
selfconsistently with electronic states by quasiclassical Eilen- 
berger theory. The selfconsistent calculation is necessary for 
quantitative estimate of physical quantities, because we have 
to determine the vortex core radius accurately. By these cal- 
culations, we investigate (1) the stable vortex lattice structure, 
(2) the local density of states (LDOS) of electrons in the vor- 
tex lattice state, and (3) the amplitude of specific heat oscilla- 
tion under field rotation. We discuss their diff'erences between 
magnetic field orientations H || node and H || antinode. In 
J;c2-);2-wave pairing, [1, 1, 0] is node-direction, and [1, 0, 0] is 
antinode-direction. 

(1) The configuration of the vortex lattice reflects the 
anisotropy of superconductivity such as J- wave pairing. 
When H II c, due to the fourfold symmetry on the Fermi sur- 
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face, triangular vortex lattice at low fields becomes the de- 
formed lattice, accompanied by a first order transition of lat- 
tice orientation change, and finally reduces to a square vortex 
lattice. When H || ^zZ?, eff'ects of anisotropic superconductiv- 
ity on the vortex lattice configuration is not enough clarified. 
Thus, in this field direction, we discuss the transition between 
two possible vortex lattice configuration, depending T or field 
orientation within the ab plane. 

(2) In the LDOS when H || c, the zero-energy electronic 
states around vortex show star-shape spatial pattern, reflect- 
ing the anisotropy of the pairing gap function. The low- 
energy states extend from the vortex core toward the node (or 
minimum gap) directions. ^^'^^^ In this paper we discuss the 
LDOS around vortices when H || ab, both for H || node and 
H II antinode. The LDOS structure mainly reflects properties 
of Q2D Fermi surface. 

(3) Since low temperature specific heat is proportional 
to zero-energy DOS N{E = 0), the oscillation of N(E = 
0) by rotating magnetic field orientation was estimated by 
Doppler shift method^^ and microscopical Eilenberger the- 
ory.^' While the Doppler shift method^^^ is a handy estimate 
to intuitively understand the oscillation, for quantitative com- 
parison with experimental data, we have to carefully evaluate 
the oscillation amplitude by quantitatively reliable theoretical 
method such as Eilenberger theory, since the oscillation am- 
plitude is in the order of a few percent of the total specific heat 
value. 

In the early stage of the specific heat experiment for 
CeCoIns,^^ the sign of the oscillation was not consistent 
with the dx2_y2-w2LYe pairing. The J;c2_^2-wave pairing was 
suggested by some experiments, such as the oscillation of 
thermal conductivity by rotating magnetic fields, flux line 
lattice transformation, ^^'^^"^^^ point-contact Andreev reflec- 
tion,^^^ and spin resonance. ^^'^^^ In order to settle the con- 
tradiction between the J;^2_^2-wave pairing and the specific 
heat oscillation, the possibility of the sign change as a func- 
tion of T was proposed.^^'^^^ Since the sign of the oscillation 
is changed at intermediate T range, the sign can be opposite 
to that expected at low T. The previous work by Eilenberger 
theory for the sign change was done by the Pesch approxi- 
mation,^^^ where the spatial dependence of the quasiclassical 
Green's function g is neglected, replacing g to the spatial av- 
erage of g. In this paper, we quantitatively calculate the am- 
plitude and sign of the specific heat oscillation by Eilenberger 
theory fully selfconsistently without using Pesch approxima- 
tion. Extending this method further, we evaluate the contri- 
bution of strong paramagnetic eff'ect and anisotropic Fermi 
surface velocity on the oscillation. The anisotropic Fermi ve- 
locity is another key factor other than the anisotropic su- 
perconducting gap, when we consider the vortex structure 
in anisotropic superconductors. The sign change of the spe- 
cific heat oscillation at low T and low H has been recently 
observed in CeCoIns, which confirms J;^;2_^2-wave pairing. 
Part of the results in this paper was reported in ref. 15 as the- 
oretical analysis for the experiment. 

It is noted that the strong paramagnetic eff'ect is impor- 
tant and indispensable when we discuss the vortex states in 
CeCoIns. The paramagnetic eff'ect comes from mismatched 
Fermi surfaces of up- and down- spin electrons due to the large 
Zeeman splitting. At higher fields in the case of strong param- 
agnetic eflfect, the upper critical field Hc2 changes to the first 



order phase transition^^'^^'^^^ and new Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) state may appear.^ ^"^^^ Even in the vor- 
tex states at low field before entering to the FFLO state, the 
strong paramagnetic eff'ects induces anomalous behavior of 
physical quantities,^^ "^^^^^ including field dependence of flux 
line lattice form factor.^^'^^^ Therefore, we consider contribu- 
tions of the strong paramagnetic eflfect in the studies of items 
(l)-(3). 

After giving our formulation of quasiclassical theory in 
the presence of the paramagnetic eff'ect in §2, we study the 
stable vortex lattice configuration in §3 and the zero-energy 
LDOS structure in §4 for two field orientations H || node and 
H II antinode, when H || ab in the (i_^2_^2-wave pairing. In 
§5, we estimate the H- and T-dependences of amplitude and 
sign of specific heat oscillation by rotating the magnetic field 
orientation within the ab-p\sine, in the presence of strong para- 
magnetic eff'ect and anisotropic Fermi velocity in addition to 
the J- wave pairing. The last section is devoted to summary 
and discussions. 

2. Formulation by Selfconsistent Quasiclassical Theory 

We calculate the spatial structure of the vortex lattice state 
by quasiclassical Eilenberger theory in the clean limit."^^"^^^ 
The quasiclassical theory is quantitatively valid when ^ » 
1/^F (^F is the Fermi wave number, and ^ is the supercon- 
ducting coherence length), which is satisfied in most of super- 
conductors in solid states. When the paramagnetic eflfects are 
discussed, we include the Zeeman term jUBB(r), where B(r) is 
the flux density of the internal field and jub is a renormalized 
Bohr magneton."^^'"^^"^^^ The quasiclassical Green's functions 
g((L>n + ifiB, k, r), f(ojn + ifiB, k, r), and /^(con + ifiB, k, r) are 
calculated in the vortex lattice state by the Eilenberger equa- 
tion 

{ojn + 1^/5 + V • (V + iA)}/ = A(f)g, 

{ojn + i^5 - V • (V - iA)} f = A'fg, (1) 

where ^ = (1 - ff^)^^^, Reg > 0, v = v/vpo, and // = 
yUB^o/^fe^c- k = (ka,h,kc) is the relative momentum of 
the Cooper pair, and r is the center-of-mass coordinate of 
the pair. We set the pairing function 0(k) = 0;^2_^2(k) = 
^/2(kl - kl)/(kl + kl) in J;^2_^2-wave pairing, and magnetic 
fields are applied to [1, 0, 0] or [1, 1, 0] directions in the crystal 
coordinate. For example, when a magnetic field is applied to 
[1, 0, 0] direction, the coordinate (x,y, z) for the vortex struc- 
ture corresponds to (b, c, a) of the crystal coordinate. In the 
case of Jjcj-wave pairing 0(k) = ^xy(^) = 2 ^f2kakb/(kl + kl), 
the results for [1,0,0] and [1,1,0] field directions are ex- 
changed. 

In our calculation, length, temperature, Fermi velocity, 
magnetic field and vector potential are, respectively, scaled 
by Rq, Tc, vf, ^0 and ^o^o- Here, Rq = fiVY^Ilnk^Tc, Bq = 
^o/2nRl with the flux quantum (po, and vp = (v^)^^ is an av- 
eraged Fermi velocity on the Fermi surface. (• • • )k indicates 
the Fermi surface average. The energy E, pair potential A and 
Matsubara frequency oji are in unit of nk^Tc. As a model of 
the Fermi surface, for simplicity, we use a Q2D Fermi surface 
with rippled cylinder- shape, and the Fermi velocity is given 
by 

vf - {va, Vb, Vc) oc (v COS 0k, V sin Gk, sin kc) (2) 
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at the Fermi surface kp = (ka^h^kc) oc 
(kj^o cos 6k, kj^o sin Ok, kc) with -n < 6k < tt and 
-71 < kc < 71?^^ To include the Fermi velocity anisotropy, we 
use V = 1 + pcosA6k. In our calculation except for last part, 
J3 = since we mainly consider the case when the Fermi 
velocity is isotropic in the ab-plsine. In our work, we set 
Vz = 0.5, so that the anisotropy ratio 



(a) type-A 



0.5, 



(3) 



as in CeCoIns. 

When magnetic fields are applied to the z axis of the vor- 
tex coordinate, the vector potential A(r) = x r + a(r) 
in the symmetric gauge, where H = (0, 0, H) is a uniform 
flux density and a(r) is related to the internal field B(r) = 
H + V X a(r). The unit cell of the vortex lattice is given by 
r = ^i(ui - U2) + S2U2 with -0.5 < Si < 0.5 (/=!, 2), 
ui = (Gx, 0, 0), U2 = (ax/ 2, Gy, 0) and GxayB = cpo. 

As for self consistent conditions, the pair potential is calcu- 
lated by 



(4) 



0<6c>„<6t>cii 



with (goNo) ^ = InT -\- 2T YjO<aj,<aj,,,^i^ - We use ^^cut = 
20kBTc. The vector potential for the internal magnetic field 
is selfconsistently determined by 

V X (V X A) = V X Mpara(r) - ^ X <^Fl^^>k ' 



0<ajn 



where we consider both the diamagnetic contribution of su- 
percurrent in the last term and the contribution of the param- 
agnetic moment Mpara(r) = (0, 0, Mpara(r)) with 



Mpara(r) = Mq 



B(r) 2T 



0<OJn 



(6) 



The normal state paramagnetic moment Mq = (jd/K)^H, k = 
Bq /jikB Tc VSttA/q and No is the DOS at the Fermi energy in the 
normal state. We set the Ginzburg-Landau parameter a: = 89 
for this typical type-II superconductor. We solve eq. (1) and 
eqs. (4)-(6) alternately, and obtain selfconsistent solutions as 
in previous works,"^^'"^^'"^^^ under a given unit cell of the vortex 
lattice. 

In Eilenberger theory, free energy is given by 

F= f drp|B(r)-Hp-yu2|5(r)p 

»^unitcell 

+|A(r)p(lnr + 2r J] cj-„') 

0<ti>„<Ct>cut 

-T Yj Wr>k,a;„))} (7) 

\a)„\<Wc„t 



with 



/(r,k,^^„) = A.^/ + A>7 

Hg - i^){7 (oJn +ifiB + y-(V + iA))/ 

+^ (w„ + + V • (V - iA)) f] (8) 
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Fig. 2. (Color online) Possible stable vortex lattice configurations type-A 
(a) and type-B (b) when H || ab. (c) and (d) show the spatial structure of 
the pair potential |A(r)| in two vortex lattice configurations a = 1.9 (a) and 
6.6 (b) for H || antinode. fi = 2, H/Hd = 0.211 and T/Tc = 0.2. Both (c) 
and (d) panels, the horizontal axis is [100] direction, and the vertical axis 
is [001] direction, and the view range is 25x25 in the Eilenberger length 
unit Rq. 



in our dimensionless unit. Using eqs. (1) and (4), we obtain 



dr{/^2|B(r-) _ n|2 . 



/unitcell 



+T Yj Re(^(A^/ + A>7)\). (9) 

The entropy in the superconducting state, given by Ss{T) = 
S„(T) - dF/dT, is obtained from eqs. (7) and (8) as 

Ss(T) 



= T-l f 

2 Juni 



dr 



^n(^c) ^ J unitcell 

(inr + 2r Z Re<A.^/ + AVY) 

0<6t>„ <6i>cut 0<a)n <Cc>cut 



A<pf + A>Y 



'+ 1 



+ 2w„Re<g-l» (10) 

\ P + I / 

0<tc>„<Ci>cut 

in our dimensionless unit. is the entropy in the normal 
state. We calculate 5's(r) numerically using selfconsistent so- 
lutions of quasiclassical Green's functions, A(r) and A(r). By 
numerical derivative of Ss(T), we obtain the specific heat C 
as 

C = T — (11) 

3. Stable Vortex Lattice Configuration When H \\ ab 

When H || ab, in the absence of paramagnetic eff'ect 
(yu = 0), the upper critical field Hc2 has field orientation de- 
pendence; Hc2 ~ 1.55 for H || antinode and Hc2 ~ 1.45 for 
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Fig. 3. (Color online) Free energy F - Fq as a function of a at T ITc = 0.1 
(a) and 0.6 (b) in the presence of strong paramagnetic effect (// = 2). 
H/Hc2 = 0.211. Fq = F(a = 1.9, H \\ node). In each figure, the upper 
(lower) line shows the case when magnetic field is applied to the node 
(antinode) direction of the J-wave pairing gap function. For each field ori- 
entation, the arrow indicates the minimum of F. 



H II node in our parameters. In the case of strong paramag- 
netic effect = 2, Hc2 is largely suppressed to Hc2 ~ 0.190 
both for H || antinode and H || node. There orientation- 
dependence of Hc2 becomes isotropic within the ab plane. 

First, we discuss the vortex lattice configuration, when 
magnetic fields are applied parallel to the ab plane, i.e., H || 
ab. In these field orientations, there are two candidates for 
the stable vortex lattice configuration. One is when one of the 
neighbor six vortices is located in the <2/7-direction, as shown 
in Fig. 2 (a), which is called type- A configuration in this pa- 
per. In the other configuration, one of the neighbor vortices is 
located in the c direction, as shown in Fig. 2 (b), called type-B 
configuration. In the isotropic case of the ^-wave pairing and 
the Fermi sphere, anisotropy ratio of the vortex lattice defined 
by 0^ = 2ay/ ^ax is a = 1 in type- A, and a = 3 in type- 
B. In our parameter where anisotropic ratio of Fermi velocity, 
~ 2, the triangular lattice is distorted by this ratio so that 
a for stable vortex lattice becomes about two times larger, 
i.e., a ~ 2 or 6. The exact value of a is not trivial because 
of the contributions of rippled cylindrical Fermi surface and 
line nodes of J- wave pairing. Thus, we have to estimate a for 
stable vortex lattice configuration. 

In Fig. 3, we plot the free energy F as a function of a for 
two magnetic field orientations dXT jT^ = 0.1 and 0.6 in the 
presence of large paramagnetic eff'ect (yu = 2). There, the lo- 
cal minimum for F is located at two configurations a ~ 1.9 
(type- A) and 6.6 (type-B). Among two local minima, vortex 
lattice configuration with lower F is stable, and the other is 
metastable. When the field is applied to the node direction of 
the J- wave pairing function, stable vortex lattice is type-A at 
T/Tc = 0.1, and changed to type-B at T/T^ = 0.6. On the 
other hand, when the field is parallel to antinode direction. 
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Fig. 4. (Color online) Temperature-dependence of the free energy diff'er- 
ence AF between two vortex lattice configurations type-A and type-B. We 
plot AF = F{a = 1.9) - F{a = 6.6), as a function of T, for two magnetic 
field orientations H || antinode (squares) and H || node (circles), ju = 2 and 
H/Hc2 = 0.211. 

Stable vortex lattice is type-B at T/Tc = 0.1, and changed to 
type-A at r/Tc = 0.6. Therefore, the vortex lattice configu- 
ration can be changed between type-A and type-B, as a first 
order transition by rotation magnetic field orientation within 
the ab plane. 

To see the temperature dependence, in Fig. 4 we show 
the free energy diff'erence AF between vortex lattice con- 
figurations type-A (a = 1.9) and type-B (a = 6.6). When 
AF > 0, type-B is stable, and when AF < 0, type-A is stable. 
From Fig. 4, we see that stable vortex lattice configuration is 
changed atT ^ 0.27^. Near T ^ 0.27^ vortex lattice configu- 
ration of type-B is stable both field orientations H || node and 
H II antinode. As seen in Fig. 3, at lower (higher) temperature 
regions stable vortex lattice configuration is type-A (type-B) 
for H II node, while the configuration of type-B (type-A) is 
stable for H || antinode. Thus, the transition of vortex lattice 
between type-A and type-B can occur upon changing temper- 
ature. 

For reference, in Fig. 5, we show the results when the para- 
magnetic eff'ect is absent (yu = 0). There, the results are qual- 
itatively unchanged from those of Figs. 3 and 4 in the strong 
paramagnetic case. The free energy as a function of a has lo- 
cal minimum at a ~ 1.9 and a ~ 6.6. At lower T, the stable 
vortex lattice configuration is type-A for H || node, and type- 
B for H II antinode. These are changed at higher T. 

The purpose of this section was that we proposed the possi- 
bility of interesting transitions in the stable vortex lattice con- 
figuration when magnetic field is applied to the ab plane. The 
vortex lattice configuration can be changed between type-A 
and type-B, depending on T, //, and the field orientation rela- 
tive to the node direction of the J-wave pairing. It is noted 
that the results of the stable vortex lattice shown in our calcu- 
lation is an example for the transition of the vortex lattice con- 
figuration. It is not sure that these results of T-dependence are 
universal. Since the free energy diff'erence AF is very small, it 
is difficult to identify the reason for changes of stable vortex 
lattice configurations. Detailed forms of Fermi surface shape 
and anisotropic pairing gap function may change the results 
for the stable vortex lattice configuration. Therefore, we have 
to carefully estimate it, using realistic Fermi surface, to com- 
pare with future experimental data. However, it is interesting 
to experimentally identify the phase diagram of the stable vor- 
tex lattice configuration under parallel field to the ab plane. 
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Fig. 5. (Color online) (a) Free energy F - as a function of a diT jT^ = 
0.6. Fq = F{a = 1.9, H \\ node). For each field orientation, the arrow in- 
dicates the minimum of F. (b) T-dependence of the free energy difference 
AF = F(a = 1.9) - F(a = 6.6), between two vortex lattice configura- 
tions type- A and type-B. Here we show the cases without paramagnetic 
effect (jd = 0) for magnetic field orientations H || antinode (squares) and 
H II node (circles). H/Hc2 ~ 0.065 with Hc2 = ^c2,antinode ~ 1-55. 



4. Electronic States in Vortex Lattice When H || aZ? 

When we calculate the electronic states, we solve eq. (1) 
with iojn E+irj. The LDOS is given by A^(r, E) = N^(r, E)+ 
Ni(v, E), where 

A^^(r, E) = No{Re{g(aJn + icr/iB, k, r)L-^„^£+,-,})k (12) 

with (T = 1 (-1) for up (down) spin component. We typically 
use T] = 0.01. The DOS is obtained by the spatial average of 
the LDOS as A^(^) = N^(E) + Ni(E) = {N(r, E))r. 

To show the electronic states in the vortex lattice for H || ab 
without the paramagnetic effect, in Fig. 6 we present zero- 
energy LDOS N(r, E = 0) for two vortex lattice configura- 
tions type-A and type-B for two magnetic field orientations 
H II node and H || antinode. There, zero-energy electronic 
states localized around vortex cores have tails extending out- 
side the core. In our calculation, since we assume an open 
Fermi surface of rippled cylinder, there are not Fermi veloc- 
ity pointing to the c-axis directions. Therefore quasiparticles 
with the a/7-direction component Fermi velocity dominantly 
contribute to the zero-energy LDOS. Therefore, N(r, E = 0) 
extends to ab direction from vortex cores, reflecting the rip- 
pled open Fermi surface. And the connection of zero-energy 
electronic states between vortices are week along the c-axis 
direction, compared to other directions. Two parallel lines of 
the LDOS tails connecting neighbor vortices are due to the in- 
terference between electrons at neighbor vortex cores. When 
zero-energy LDOS are connected between the vortices, the 
LDOS is suppressed just on the straight line directly connect- 
ing the vortex centers, resulting in two parallel tails of the 
zero-energy LDOS.^^^ As for the difl'erence between H || node 




Fig. 6. (Color online) Spatial structure of zero-energy LDOS N(E = 0, r) 
for // = at H/Hc2 = 0.065 in the vortex lattice under magnetic fields 
parallel to the ab plane. Upper two panels (a) and (b) are for a magnetic 
field along the antinode direction of the J-wave paring gap function. The 
vertical axes are [001] directions and the horizontal axes are [100] direc- 
tions. Lower two panels (c) and (d) are for the node-direction. The vertical 
axes are [001] directions and the horizontal axes are [110] directions. Left- 
side two panels (a) and (c) are for the vortex lattice configuration of type-A 
(a = 1.9). Right-side two panels (b) and (d) are for type-B configuration 
(a = 6.6). The view ranges of all four panels are 15x15 in the Eilenberger 
length unit Rq. 



and H || antinode, zero-energy states around vortex core are 
broadly extended towards outside when H || node, because su- 
perconducting gap has node when quasi-particles propagate to 
the direction perpendicular to the magnetic field orientation. 
Zero energy LDOS in the presence of strong paramagnetic ef- 
fect (yu = 2) is presented in Fig. 7. Due to the Zeeman splitting 
by strong paramagnetic eff'ect, the original zero-energy bound 
state for yu = is lifted to finite energy.^^^ Thus, the tails ex- 
tending from vortex cores are smeared in zero-energy LDOS 
of Fig. 7. 

5. Field-Angle Dependence of Specific Heat 

Using self consistent results of A(r), A(r) and quasiclassi- 
cal Green's functions, we calculate the entropy Ss(T) by eq. 
(10). In Fig. 8, we show T-dependence of Ss(T). From the 
T-dependence, we numerically obtain T-dependence of the 
specific heat by eq. (11), which is presented in Fig. 9. There, 
C oc at low fields because of line nodes in the J-wave 
pairing. As shown in Fig. 9(a) when paramagnetic eff'ect is 
negligible (pi = 0), with increasing H, C reduces to T-linear 
behavior due to low energy excitations around vortices, and 
approaches the line for the normal states. Since Tc decreases 
with increasing H, the jump of C at Tc becomes smaller at 
high fields. When paramagnetic eff'ect is strong (yu = 2) as 
shown in Fig. 9(b), C does not deviate from the low field 
curve even at H/Hc2 = 0.42, because H(= 0.08) is still small. 
It changes to the normal state by the first order transition at 

Hc2- 
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Fig. 7. (Color online) The same as in Fig. 6, but in the presence of strong 
paramagnetic effect (ji = 2) and H/Hc2 = 0.211. The view ranges of all 
four panels are 25x25 in the Eilenberger length unit Rq. 




Fig. 9. (Color online) T-dependence of Specific heat C at H/Hc2 =0.065, 
0.194, 0.323, and 0.452 from bottom to top for // = (a), and at 
H/Hc2 =0.105, 0.211, 0.316, and 0.421 from bottom to top for // = 2 
(b). Dotted line indicates C(T) in the normal state. 





Fig. 8. (Color online) T-dependence of entropy S s(T)/Sn(Tc) at 
H/Hc2 =0.105, 0.211, 0.316, and 0.421 from bottom to top. ju = 2. 



Next, we discuss the oscillation of the specific heat in the 
form 



C = Ch(1-A4 cos 4^), 



(13) 



when magnetic field orientation is rotated within aZ?-plane. 
^ is a relative angle between field orientation and antinode- 
direction of the J-wave gap function. The oscillation factor 
A4 is calculated by 

. ^node ~ ^antinode 



X 100 [%] 



(14) 



^node "I" ^antinode 

where Cnode (Qntinode) IS the specific heat for field orientation 
along the node (antinode) direction, i.e., 6 = 45° (0°) in the 
(i;c2 -3/2 -wave pairing. 

To evaluate how A4 depends on the vortex lattice configura- 
tions type- A and type-B, in Fig. 10 we present T-dependence 
of A4 for some choices (a^node, Ofantinode) = (1.9,1.9), (6.6,6.6), 

(1.9,6.6), (6.6,1.9). <^node (^^antinode) is ValuC of a fOY VOrtCX 



Fig. 10. (Color online) T-dependence of fourfold oscillation part A4 in spe- 
cific heat under magnetic field rotation at /u = 2 and H/Hc2 = 0.211. We 
show A4 for some vortex lattice configurations (o^node. Qfantinode) =(6.6,6.6), 
(1.9,6.6), (6.6,1.9), (1.9,1.9) from bottom to top. 



lattice configuration when field orientation is along the node 
(antinode) direction. From the figure, we see that the choice 
of the stable vortex lattice configuration does not seriously af- 
fect on the behavior of A4. In any cases, on lowering T from 
high temperature, A4 increases with sign change from nega- 
tive to positive, and after peak near T ~ O.STc, A4 rapidly 
decreases with sign change to negative. Therefore, we fix the 
vortex lattice configuration as type- A with a = 2, hereafter. 
We note that C4 for a = 2 does not seriously change from that 
for a = 1.9. 

Figure 1 1(a) shows T-dependence of A4 for some H. There, 
for lower //, sign change of A4 occurs at T ~ OAST^. How- 
ever, at higher H, A4 remains positive and increases on low- 
ering T. To show the H- and T-dependence, we show con- 
tour plot of A4(T,H) in Fig. 11(b). There solid line indicate 
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Fig. 11. (Color online) Fourfold oscillation part A4 of specific heat in the 
case of strong paramagnetic eff'ect = 2. (a) T-dependence of oscillation 
part A4 at HIHci =0.105, 0.211, 0.316, and 0.421. (b) (7, ^)-dependence 
of A4. The solid curve indicate position of A4 = where sign of the oscil- 
lation changes. ~ 0.19.yS = 0. 



the sign change of A/^{T, H). Inside region from the solid line 
at lower T and lower H has negative A4. Other higher T re- 
gion has positive A4. We defined the magnetic field where the 
sign change occurs at low temperature as ifs, and the temper- 
ature of the sign change as T^. In Fig. 1 1(b) in the presence of 
strong paramagnetic eff'ect, T^IT^ -0.1 and H^IHd ~ 0.35 
[H, ~ 0.067]. 

To evaluate the contribution of Pauli-paramagnetic eff'ect, 
in Fig. 12 we show A/^iT, H) in the case without paramagnetic 
eff'ect (yu = 0). There, sign change of A4 similarly occurs at 
low T and //, while amplitude of A4 is enhanced compared 
with that for yu = 2 in Fig. 11. This case (yu = 0) corresponds 
to the previous work in ref. 27. Our results are qualitatively 
consistent to it, while we perform selfconsistent calculation 
without using Pesch approximation. Quantitatively, Hs/Hc2 ~ 
0.35 in our calculation becomes smaller than H^IHd ~ 0.5 in 
ref. 27. T^IT^ is almost same in both calculations. The reason 
of the sign change is because of the £^-dependence of DOS 
N{E), as discussed in ref. 27. Compared with N{E) for H || 
node, N{E) for H || antinode is smaller at low but larger at 
higher E. 

Our new approach is to evaluate the contribution of para- 
magnetic eff'ect, by comparing results in Figs. 11 and 12. By 
the paramagnetic eff'ect, the upper critical field is largely sup- 
pressed from Hc2 ~ 1.55 (when yu = 0) to //c2 ~ 0.19 (when 
yu = 2), and the sign change field is also largely suppressed 
from H, ~ 0.54 to H, ~ 0.067. Thus, in the scale HjHa, 
the normalized sign-change field H^IH^i keeps similar value, 
which is not seriously changed by the paramagnetic eff'ect. 
We see also the diff'erence at // > H^/m Fig. 11, A4 mono- 
tonically increase on low temperature, but in Fig. 12, A4 de- 
creases at low temperature after increase at high temperature. 



Fig. 12. (Color online) The same as Fig. 11, but in the case when the para- 
magnetic eff'ect is absent {ji = 0). In (a), H/Hc2 =0.065, 0.194, 0.323, and 

0.452. H,2 ^ //c2,antinode ~ 1-55. = 0. 



Thus, solid line of //§ almost horizontal at low T in Fig. 1 1(b), 
but the line of //§ increases on lowering T in Fig. 12(b). 

Lastly, we also evaluate the contribution of Fermi veloc- 
ity anisotropy on A4, considering finite in the definition of 
Fermi surface in eq. (2). Figure 13 shows A4(T, H)for JS = 0.5 
and jd = 2. This positive JS shifts smaller to ifs/^c2 ~ 0.3 
foryS = 0.5, from HJHci = 0.35 foryS = [Fig. 11]. On the 
other hand, negative j3 makes //§ larger. Therefore, anisotropic 
Fermi velocity aff'ect on the sign-change field H^. However, 
from Figs. 1 1 - 13, we see that the sign-change temperature 
of A4 is rather independent from the Fermi surface anisotropy 
and paramagnetic eff'ect. 

Recently the sign-change of A4 in specific heat oscillation 
was observe in CeCoIns.^^^ There, A4 shows the sign change 
at low T and low H. This is qualitatively consistent to our re- 
sults of calculations, and supports the d^i.yi -waye pairing for 
superconductivity in CeCoIns. Compared to the experimental 
data in ref. 15, amplitude of A4 is smaller in our calculation. 
One of the reason is that in analysis of experimental data def- 
inition of A4 is C = Co + C//(l - A4COs4^), while in our 
calculation Co = in the clean limit. While T^/Tc ~ 0.1 
both for experiment and theory, H^/Hc2 ~ 0.1 in experi- 
mental data is smaller than the theoretical calculation. We 
showed that //s///c2 can be changed by the Fermi velocity 
anisotropy, based on a simple yS-model. Therefore, as a pos- 
sibility, by theoretical estimate using realistic Fermi surface 
structure may improve quantitative accordance of with ex- 
perimental data. 

6. Summary and Discussions 

We investigated the vortex state in (i^2_^2-wave pairing 
when magnetic field H is applied parallel to ab plane, based 
on quantitative calculation by selfconsistent Eilenberger the- 
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Fig. 13. (Color online) The same as Fig. 11, but in the case when we in- 
clude the contribution from anisotropic Fermi velocity hy (3 = 0.5 with 
strong paramagnetic effect // = 2. In (a), /////c2 =0.109, 0.217, 0.326, and 
0.435. //c2 ~ 0.184. 

ory. Evaluating (1) stable vortex lattice structure, (2) the zero- 
energy LDOS in the vortex lattice state, and (3) temperature 
dependence of specific heat, we estimate the diff'erences of 
vortex states for H || node and for H || antinode. 

(1) The transition of two possible vortex lattice configura- 
tions (type- A and type-B in Fig. 2) can occur as a function 
of temperature and magnetic field. There is a possibility that 
vortex lattice configuration for H || antinode is diff'erent from 
that for H || node. This indicate that by rotation of magnetic 
field orientation within ab plane, first order transition occurs 
between two vortex lattice configuration. 

(2) The spatial structure of zero-energy LDOS reflects Q2D 
Fermi surface structure. The tails of zero-energy LDOS pre- 
fer connecting with those of neighbor vortices in the ah di- 
rection compared to those in c direction. The dependence on 
the relative angle between the node-direction and magnetic 
field orientation gives minor contribution to broad extension 
of quasiparticles around vortex cores. 

(3) We evaluated magnetic field and temperature depen- 
dence of the amplitude and sign for specific heat oscillation 
by rotation of magnetic field orientation. The sign of the os- 
cillation changes at low field and low temperature region, 
where the sign is consistent to the oscillation of evaluation 
by zero-energy DOS. This sign-change behavior was recently 
observed in CeCoIns.^^^ Our selfconsistent calculation with- 
out Pesch approximation gives qualitatively consistent results 
to those in previous work by Pesch approximation when the 
paramagnetic eff'ect is absent.^^^ We extend this calculation to 
the case of strong paramagnetic eff'ect such as CeCoIns, and 
show the sign-change behavior of the specific heat oscillation 
is similar to that without paramagnetic eff'ect, if the magnetic 
field is scaled by suppressed //c2. We also evaluate the con- 



tribution of anisotropic Fermi velocity, which shifts the sign- 
change field. 

The experiments to observe the dependence on magnetic 
field orientation are useful method to identify the position of 
the node in the superconducting gap on the Fermi surface. In 
addition to specific heat oscillation by rotation of magnetic 
field orientation, we examine the dependence on the relative 
angle between the node and field orientation for stable vortex 
lattice configuration and spatial structure of LDOS, with and 
without paramagnetic eff'ect. The theoretical calculations give 
valuable information to be compared with details of experi- 
mental data. 
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